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Abstract 

The equation of state for quark matter is derived for a nonlocal, chiral quark model 
within the mean field approximation. We investigate the effects of a variation of 
the formfactors of the interaction on the phase diagram of quark matter. Special 
emphasis is on the occurrence of a diquark condensate which signals a phase tran- 
sition to color superconductivity and its effects on the equation of state under the 
condition of (3- equilibrium and charge neutrality. We calculate the quark star con- 
figurations by solving the Tolman- Oppenheimer- Volkoff equations and obtain for 
the transition from a hot, normal quark matter core of a protoneutron star to a cool 
diquark condensed one a release of binding energy of the order of AMc 2 ~ 10 53 erg. 
We find that this energy could not serve as an engine for explosive phenomena since 
the phase transition is not first order. Contrary to naive expectations the mass de- 
fect increases when for a given temperature we neglect the possibility of diquark 
condensation. 
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1 Introduction 



Color superconductivity (CS) in quark matter [1] is one interesting aspect of 
the physics of compact star interiors [2]. Since calculations of the energy gap 
of quark pairing predict a value A ~ 100 MeV and corresponding critical tem- 
peratures for the phase transition to the superconducting state are expected 
to follow the BCS relation T c = 0.57 A, the question arises whether diquark 
condensation can lead to remarkable effects on the structure and evolution 
of compact objects. If positively answered, CS of quark matter could provide 
signatures for the detection of a deconfined phase in the interior of compact 
objects (pulsars, low-mass X-ray binaries) along the lines of previously sug- 
gested strategies [3,4,5,6,7]. 

Diquarks gained popularity in astrophysics about a decade ago, when they 
were suggested to influence the supernova collapse and "bounce-off" [8,9,10,11,12], 
and to enhance the neutrino cooling of quark-stars. The latter effect is now 
subject to much research within improved scenarios Refs. [13,14,15]. Possible 
consequences for the magnetic field and glitch-type phenomena are discussed 



Hong, Hsu and Sannino have conjectured [19] that the release of binding energy 
due to Cooper pairing of quarks in the course of protoneutron star evolution 
could provide an explanation for the unknown source of energy in supernovae, 
hypernovae or gamma-ray bursts, see also [20] . Their estimate of energy release 
was based on the simple formula 



which did not take into account the change in the gravitational binding energy 
due to the change in the structure of the stars quark core. In the present work 
we will reinvestigate the question of a possible binding energy release due to 
a CS transition by taking into account changes in the equation of state (EoS) 
and the configuration of the quark star selfconsistently. 

As a first step in this direction we will discuss here the two flavor color super- 
conducting (2SC) quark matter phase which occurs at lower baryon densities 
than the color-flavor-locking (CFL) one, see [21,22]. Studies of three-flavor 
quark models have revealed a very rich phase structure (see [22] and refer- 
ences therein). However, for applications to compact stars the omission of the 
strange quark flavor within the class of nonlocal chiral quark models consid- 
ered here may be justified by the fact that central chemical potentials in stable 
star configurations do barely reach the threshold value at which the mass gap 
for strange quarks breaks down and they appear in the system [23] . Therefore 



in [16,17,18]. 
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we will not discuss here first applications to calculate compact star configu- 
rations with color superconducting quark matter phases that have employed 
non-dynamical quark models [24,25]. Their results depend on the somewhat 
arbitrary choice of the strange quark mass m s , the diquark gap A, and the 
bag constant. Dynamical calculations like the one we are going to present here 
fix their model parameters with the constraint to describe observable hadron 
properties. They include a selfconsistent solution of the chiral and diquark 
gap equations which shows that the bulk of a compact star interior does not 
contain strange quark matter [23,26]. 

We will investigate the influence of the formfactor of the interaction on the 
phase diagram and the EoS of dense quark matter under the conditions of 
charge neutrality and isospin asymmetry due to /9-equilibrium relevant for 
compact stars. 

Finally we consider the question whether the effect of diquark condensation 
which occurs in the earlier stages of the compact star evolution (t ~ 100 s) 
[15,27,28] at temperatures T ~ T c ~ 20 — 50 MeV can be considered as an 
engine for exposive astrophysical phenomena like supernova explosions due to 
the release of a binding energy of about 10 52 10 53 erg, as has been suggested 
before [19,20]. 



2 Thermodynamic potential for asymmetric 2SC quark matter 

We consider a nonlocal chiral quark model described by the effective action 

s$, v>]=E W - ™)<l> + S *S, *i>] + s m [A # (2) 

P 

S q q[*P,^] = G 1 J2 [{i>g{p)^){$g{p')i>) + (ipil5Tg(p)^)(*pii5Tg(p')ip)] , (3) 
p,p' 

S„$, ifl = G 2 Y,(^ 2 X 2 Cg(p)^ T )^ T Ct l5 r 2 X 2 g(pm, (4) 
p,p' 

which generalizes the approach of Ref. [29] by including the scalar diquark 
pairing interaction channel with a coupling strength G 2 . Here, the matrix 
C = 27072 is the charge conjugation matrix operator for fermions, the matri- 
ces t 2 and A 2 are Pauli and Gell-Mann matrices acting in, respectively, the 
SU(2) flavor and SU(3) color spaces; rh is the diagonal current quark mass ma- 
trix m = diag(m u , m d ). The function g{p) is the formfactor of the separable 
four-fermion interaction. In the present work it is assumed to be instanta- 
neous, so that it depends only on the modulus of the three-momentum. After 
bosonization and integrating out the quark fields, the mean field approxima- 
tion introduces the following order parameters: the diquark gap A, which can 
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be seen as the gain in energy due to diquark condensation, and the mass 
gaps (f) u , <pd, which indicate chiral symmetry breaking. The grand canonical 
thermodynamic potential per volume can be obtained as 



/d?p 1 / 1 ~ \ 

(27r)3 2 Trln \ f S l<yiuJn ^J 



(2tt) 3 2 



8Gi 



4 Go 



(5) 



where oo n = (2n+l)7iT are the Matsubara frequencies for fermions. The chem- 
ical potential matrix ft = diag(// u , /id) as well as the dynamical quark mass 
matrix M = diag(m u + <p u g(p),md + 4>dg{p)) are included into the definition 
of the inverse quark propagator 



' i> - M - /t7o A7 5 T 2 X 2 g(p) 
K -A*-f 5 T 2 X 2 g{p) i> - M + /t7o 



(6) 



One can combine the chemical potentials fi u , fid of u and d quarks by intro- 
ducing fi q = (fi u + /id)/2 and fii = (fi u — /id)/2 as the Lagrange multipliers 
related to, respectively, the quark number density n q and the isospin asym- 
metry rij. One can easily check the relations n q = n u + rid = —dQ/dfi q , 
ni = n u — n d = — dVt/diii. Following Ref. [30] we use in the following that 
the chiral gaps of both flavors are degenerate (p u = <Pd = <P (and therefore also 
the masses m u = rrid = m). The trace operations in the bispinor, flavor and 
color spaces as well as the Matsubara summations can be performed with the 
following result for the thermodynamic potential of asymmetric quark matter 
[31,32,33] 



6 2 A 2 
Q g (<j>, A; n q , (j, i, T) — -t-^t + — ^ 2 
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The number of color degrees of freedom in our model is fixed to be N c = 3. 
Only the coupling constant G\ is considered to be a free parameter of the 
model, while G 2 could be chosen according to Ref. [34] as G 2 = G\ N c / (2N C — 
2). We discuss later the special choices G 2 = 0.75 G 1 and G 2 = G 1 . 



The dispersion relation — y q 2 + (m + <fi g(q)) 2 holds for unpaired quarks 
also in color-symmetry broken 2SC phase. The quarks with paired colors have 
a modified dispersion relation 



E ± = (E 4> ± N ) 



1 A 2 9 2 (q) 



For nonvanishing A in Eq. (7) the color symmetry is broken. Two of the three 
quark color degrees of freedom are coupled to bosonic Cooper pairs in the color 
antitriplet state which can form a Bose condensate. The dispersion relation 
for quarks of these colors contains the diquark gap A and becomes 

^ = (^± W ) J + AV(!) • (8) 



For quarks of the remaining third color there is no pairing gap in their disper- 
sion relation which reads 



E+ f = Jq' + imf + ^f g(q)) 2 . (9) 



In thermal equilibrium characterized by a fixed set of thermodynamic variables 
{a*u, l^d, T}, the thermodynamic potential Q q shall attain a global minimum 
with respect to all its order parameters. In the corresponding state the gap 
equations for the order parameters are fulfilled. However, the opposite is not 
always true. Satisfying the gap equations is a necessary, but not a sufficient 
condition for the thermodynamic equilibrium state. 

In Eq. (7), we subtract the term f2™ c in order to define the thermodynamic 
potential such that the vacuum pressure vanishes, 

P q (0, 0, 0) = . (10) 



On the mean-field level the present approach is analogous to a flavor- asymmet- 
ric generalization of the instanton-motivated model of Ref. [34] for two-flavor 
color superconductivity (2SC). 

From BCS theory it is known, that in order to form Cooper pairs at T = in a 
dense Fermi system, the difference in the chemical potentials of the Fermions 
to be paired should not exceed the size of the gap. As previous calculations 
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within this type of models have shown [34] , there is a critical chemical potential 
for the occurrence of quark matter /// > 300 MeV and values of the gap in 
the region A < 150 MeV have been found. Therefore it is natural to consider 
the problem of the color superconducting (2SC) phase with the assumption, 
that quark matter is symmetric or very close to being symmetric (/i u ~ fid)- 

The number densities derived from thermodynamic potential can be expressed 
in an integral form as 



00 2 J 

n q ( N , m, T) = 2 f {f(E^ -ti q - tn, T) + f(E^ - N + m, T) 
o 

— f(E<t> + Vq- 1*1, T ) - f( E <t> + Vq + Vl, T) 

-2 dE + /dfx q [f(E + - m, T) + f(E + + m, T) - 1] 

-2 dE./dfi q [f(E. - m,T) + /(£_ + m,T) - 1]} (11) 



ni(fi q , pi, T) = 2 J {f{E 4> - fi q - m, T) - f{E 4> - fx q + m, T) 
o 

+f(E 4> + fj, q - m, T) - f{E^ + fig + m, T) 
+2 \f(E + -ii I ,T)-f(E+ + ii I ,T) 

+ f(E_ - m, T) - f(E_ + fij, T)]} . (12) 

Here, f(E,T) = [1 + exp(E /T)] _1 is the Fermi distribution function and the 
derivatives are given by dE±/d/i q = (fj, q ± E^)/E±. 

Using the entropy density s q = —dQ q /dT, the pressure P q = — Q q and the 
thermodynamic identity P q + e q = s q T + fi q n q + fiifij, the EoS for the energy 
density as a function of temperature and chemical potentials becomes 



00 2 J 

/q dq 



+f(E 4> + fi q - m, T) + f{Et + iiq + T) - 2] 

+2 (E + - n q dE+/dn q ) [f(E + - m, T) + f(E + + N , T) - 1] 

+2 (E. - n q dE_/d N ) [/(£_ - m,T) + f{E. + m,T) - 1]} 

In the next sections we will apply this model for the construction of the EoS 
for stellar matter in /3-equilibrium. 
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3 EoS for 2SC quark matter in /^-equilibrium 



We consider stellar matter in the quark core of compact stars consisting of 
electrons in chemical equilibrium with u and d quarks. Hence fid = fi u + fi e , 
where fi e is the electron chemical potential. The thermodynamic potential of 
such matter is 

A; fx q , fj,j, /i e , T) = n q (<f>, A; fx q , fxj, T) + Q e (fx e , T), (14) 



with 



n e (fx e ,T) = 



12tt 2 



1 7 

r* T ~ i8o 



(15) 



for ultrarelativistic electrons. The number density of electrons is then n e = 
—dVL/dn e . The condition of local charge neutrality 

2 1 

-n u - -n d - n e = 0, (16) 



in the new notation reads 

n q + 3nj — 6n e = 0. (17) 



In our model we assume the neutrino chemical potential to be zero, which 
corresponds to a scenario where the neutrinos (and/or antineutrinos) leave 
the star as soon as they are created. The limits of such an approximation 
could be based on the calculations for the case of trapping of neutrinos in the 
hot stars interior [21]. For possible consequences of antineutrino trapping on 
hot quark star evolution, see [35]. 

In [36,32] it is shown that at given fx q the diquark gap is independent of the 
isospin chemical potential for \fii{fi q )\ < fii c (fi q ), otherwise vanishes. Increase 
of isospin asymmetry forces the system to undergo a first order phase transition 
by tunneling through a barrier in the thermodynamic potential (14). In order 
to describe this transition we have to choose the apropriate set of Lagrange 
multipliers. Under the constraints of (5— equilibrium and charge neutrality, the 
free enthalpy density can be written in the form 

G = fi u n u + jj, d n d + fi e n e — fi B n B , (18) 



where fi B = fiu + ^d is the chemical potential conjugate to the baryon number 
density n B = (n u + n<j)/3 and remains as an independent thermodynamical 
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variable under the given constraints. Note that ^ ^ A*g/3 because of the 
isospin violation in weak interactions, which determine the equilibrium condi- 
tions for compact stars. This is in contrast to equilibrium conditions for dense 
matter in heavy-ion collisions. The corresponding set of equations of state is 
obtained from the thermodynamic potential (14) in the form 



P(T,» B ) = -Sl(T,» B ) , (19) 

n B {T, im b ) = -dtt{T, fi B )/dfi B , (20) 

s(T,fi B ) = dn(T,fi B )/dT , (21) 

e(T, n B ) = T s(T, fx B ) - P(T, fi B ) + fi B n B (T, /x B ) , (22) 



where due to the above constraints the dependence on the set of thermody- 
namic potentials (fj, q , /ij, /i e ) is reduced to the dependence on /i B only. 

In concluding this section we would like to comment on the order of the CS 
phase transition within above approach. The phase transition to the 2SC phase 
coming from the hot normal state of matter by cooling the system occurs as a 
second order phase transition without a jump of concentrations, whereas the 
phase transition from an isospin asymmetric state to a superconducting one 
by decreasing the asymmetry is of first order and has a jump in the densities 
(in particular in the electron fraction). 



4 Results of model calculations 

4-1 EoS of quark matter in 2SC phase for finite temperature 

The gaps have been calculated by minimizing the thermodynamic potential 
Eq. (14), in the space of order parameters (0, A) and the results are shown 
in Figs. 1 and 2. 

In order to compare the effects of the quark interaction formfactors on the EoS, 
we study the Gaussian (G), Lorentzian (L) and cutoff (NJL) type formfactors 
defined as 

g L (q) = [l + (q/A L ) 2a }-\ a > 1 , (23) 
^ G (g) = exp(-g 2 /A^) , (24) 
3NjL(g) = #(l-g/ANjL) • (25) 

The Lorentzian (L) interpolates between a soft (Gaussian-type, a ~ 2) and a 
hard cutoff (NJL, a > 30) depending on the value of the parameter a. 
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We will employ parametrisations of the nonlocal quark model which reproduce 
pion properties: mass = 140 MeV, decay constant f n = 93 MeV and 
which have the same quark mass gap <f>(T = 0,/x = 0) = 330 MeV in the 
vacuum. The results for the parameterization are taken from [29], see Tab. I. 
G 2 remains a free parameter of the approach. For G 2 = Gi, the 2SC phase 
with nonvanishing diquark condensate is obtained whereas for G 2 = 0.75 G\ 
no 2SC phase is possible under neutron star conditions. 

Figure 1 displays the T = solutions of the chiral gap (f> and diquark gap A for 
different formfactors in f3— equilibrium for G 2 — G±. For the densities relevant 
for stable star configurations, fi q < 450 MeV, the critcal chemical potential 
\x c for the chiral transition and for the onset of diquark condensatuion does 
depend on the type of the formfactor. The value of the diquark gap reaches 
A ~ 200 MeV in the range of densities relevant for compact stars. It is rather 
insensitive to changes of the formfactor of the quark interaction. 

Figure 2 shows the solutions for the chiral gap <fi and the diquark gap A 
for the Gaussian model formfactor at different temperatures T = 0, 40, 60 
MeV. We compare results in the chiral limit (mo = 0) with those for finite 
current quark mass m = 2.41 MeV and observe that the diquark gap is 
not sensitive to the presence of the current quark mass, which holds for all 
formfactors in Eqs. (23)- (25). However, the choice of the formfactor influences 
the critical values of the phase transition as displayed in the quark matter 
phase diagram (/i q — T plane) of Fig. 3, see also Fig. 1. A softer formfactor 
in momentum space gives lower critical values for T c and ji c at the borders of 
chiral symmetry restoration and diquark condensation. The numerical results 
for the critical chemical potentials at T = are summarised in Tab. 1 for the 
different formfactors. 

The inset of Fig. 3 shows that the generalisation of the BCS relation T c ~ 
0.57 A(T = 0,/ig) g(fj, g ), between the critical temperature T c of the super- 
conducting phase transition and the pairing gap A at T = is satisfactorily 
fulfilled in the domain of the phase diagram relevant for compact stars. 

Figure 4 shows the resulting equation of state (EoS) for quark matter, as 
pressure P versus baryochemical potential hb for the isotherms T = 0,40,60 
MeV and in /^-equilibrium. The Gaussian formfactor has been used. The Fig- 
ure displays the defect of nonconfining quark models: at nonzero temperature 
quasifree quarks contribute to the pressure at low chemical potentials /is (cor- 
responding to energy densities e < 300 MeV/fm 3 ). The dot-dashed curve in 
Fig. 4 shows that the presence of a diquark condensate stabilizes the quark 
matter system and entails a hardening of the EoS. 

In Figure 5 we show the matter distribution inside the quark star and the ac- 
cumulated mass for quark stars with the same central baryon density tib(0) = 
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5 no for different temperatures. We have cut the density distribution at finite 
temperatures for n# < 0.5 no due to the following reason. At nonzero tempera- 
tures the mass gap decreases as a function of the chemical potential already in 
the phase with broken chiral symmetry. Hence the model here gives unphysical 
low-density excitations of quasi-free quarks. A systematic improvement of this 
situation should be obtained by including the phase transition construction to 
hadronic matter. However, in the present work we circumvent the confinement 
problem by considering the quark matter phase only for densities % > 0.5 n , 
where n is nuclear saturation density. 

4-2 Configurations of hot quark stars 

Here we compare configurations with and without CS, in order to investigate 
the effect of diquark condensation on the total energy of a quark star and to 
decide whether the corresponding phase transition could serve as an engine for 
explosive phenomena like supernova? and gamma ray bursts. The CS transition 
can occur during the cooling evolution of a hot protoneutron star [15,27]. We 
suppose that the approximation of an isothermal temperature distribution [40] 
in the star interior is sufficient for the estimate of the mass defect 2 . In the 
present paper we consider selfbound configurations of pure quark matter as 
models for quark cores in hybrid stars which should be the general case to be 
discussed in a separate paper, see e.g. Ref. [41] for T = configurations. The 
extrapolation from pure quark star configurations to hybrid ones, however, is 
strongly model dependent and cannot be done without a calculation. Quark 
matter effects are expected to be still valid in hybrid stars when rescaled with 
a monotonously rising function of the ratio of the quark core volume to the 
quark star volume. 

The spherically symmetric, static star configurations are defined by the well 
known Tolman-Oppenheimer-Volkoff equations [38] for the mechanical equi- 
librium of self-gravitating matter (see also [37,39]) 



Here e(r) is the energy density and P(r) the pressure at the distance r from 
the center of the star. The mass enclosed in a sphere with radius r is defined 



alP{r) 



dr 



[e(r) + P(r)} [m(r) + 47rr 3 P(r)] 
r[r — 2m(r)] 



(26) 



by 



r 




(27) 



o 



2 For a discussion of isentropic hot quark star configurations, see [37]. 
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These equations are solved with the boundary condition of a given central 
baryon number density n#(0) defining a configuration. In Fig. 6 we show 
the results for the dependence of the masses on the central energy density 
as well as on the radius of the configurations for different temperatures and 
two formfactor models: a Gaussian and a Lorentzain one. Stable configurations 
correspond to the rising branch on the mass-central density plane. In Tab. 2 we 
give for fixed temperatures T = 0, 40, 60 MeV the values of maximal possible 
mass, corresponding radius, central baryon number density and energy density 
for a Gaussian and a Lorentzian formfactor, respectively. For comparison, we 
have included T = configurations for a small coupling constant G 2 = 0.75 G 1: 
for which the diquark condensate does not occur. The comparison shows that 
the maximum mass decreases with increasing temperature due to a softening 
of the EoS. Also the disappearance of the diquark condensate (in our case due 
to lowering of the coupling G 2 ) softens the EoS and lowers the maximum mass. 
The Gaussian formfactor has the lowest critical baryochemical potential when 
compared with the Lorentzian and the NJL model. Therefore the Gaussian 
quark star configuration can extend to larger radii where the density is lower 
and can thus carry more mass than the Lorentzian model with the same central 
baryon density, see Fig. 6. 

Since in our calculations we do not consider the existence of a possible hadronic 
shell of the star where the transition to the phase with broken chiral symmetry 
and confinement occurs, we introduce a phenomenological cutoff in the baryon 
number density at the value n B = 0.5 n , see Fig 5. With this condition we 
define the radius of the quark core R. The total mass M of the object is defined 
by the sphere with that radius M = m(R). It is easy to see from that Figure 
that temperature plays no role for the mechanical equilibrium of the quark 
core close to the center of the configuration due to its high density. It has an 
influence, however, on the peripheral region and changes the total mass for a 
given central density, see Figs. 6. We have checked numerically, that a change 
of the somewhat arbitrary cutoff density does not change the qualitative result 
for the energy release reported below. 

For the discussion of evolutionary variations of the internal structure during 
the cooling it is important to take into account the conservation of the total 
baryon number of the star. Therefore, one has to compare configurations with 
different central densities, but with the same baryon number Nb, see Figs. 7, 
8. 

Here we will discuss two scenarios for the protoneutron star cooling which we 
denote by A and B, where A stands for cooling of a star configuration with SC 
whereas B is a scenario without SC. The initial states for both scenarios are 
chosen to have the same mass Mi(A) = Mi(B) for a given initial temperature 
of T = 40 MeV. The final states at T = 0, however, have different masses 
Mf(A) 7^ Mf(B) while the total baryon number is conserved in the cooling 



11 



evolution. The resulting mass differences are AM (A) = 0.06 M , AM(B) = 
0.09 M and AM (A) = 0.05 M , AM{B) = 0.07 M for the Gaussian and 
Lorentzian models, respectively. 

We make the observation that, contrary to naive expectations of the mass 
defect as a result of the binding energy from Cooper pairing [19], AM c 2 ~ 
(A/ /i) 2 M c 2 oc 10 53 erg, the mass defect in the cooling evolution with diquark 
condensation (A) is about 40% smaller than in the scenario without SC (B), 
regardless of the choice of formfactor. The difference amounts to ~ O.O2M 
for initial configurations with masses of about 1.4 — 1.5 M . Therefore, we 
conclude that general relativity effects play an important role in the estimate 
of the mass defect since the effects of diquark condensation on the mass of 
configurations via the change in the equation of state overcompensate effects 
due to the occurence of the diquark gap in the one-particle energies. 



5 Conclusion 

We have investigated the influence of diquark condensation on the thermo- 
dynamics of quark matter under the conditions of /3-equilibrium and charge 
neutrality relevant for the discussion of compact stars. The EoS has been 
derived for a nonlocal chiral quark model in the mean field approximation, 
and the influence of different formfactors of the nonlocal, separable interac- 
tion (Gaussian, Lorentzian, NJL) has been studied. The model parameters are 
chosen such that the same set of hadronic vacuum observables is described. 
We have shown that the critical temperatures and chemical potentials for the 
onset of the chiral and the superconducting phase transition are the lower the 
smoother the momentum dependence of the interaction formfactor is. 

The phase transition to color superconducting quark matter from the lower 
density regions at small temperatures (T < 5 4- 10 MeV) is of first order, 
while the melting of the diquark condensate and the corresponding transition 
to normal quark matter at high temperatures is of second order. The flavor 
asymmetry in neutral quark matter increases when the electron fraction drops. 

When the diquark coupling constant is large enough (G2 ~ G±) then the 
asymmetry which is present in the system due to the /9-equilibrium condition 
does not destroy the diquark condensation. The electron fraction under these 
conditions does not exceed n e /n total ~ 0.01. The masses of the quark core 
configurations could be up to 1.8 M for the Gaussian and up to 1.6 M 
for the Lorentzian formfactor models. The quark core radii with Lorentzian 
formfactors are up to 9 km, while for Gaussian models could extend up to 11 
km. 
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Diquark condensation makes the EoS harder, which leads to an increase in the 
maximum mass of the quark star configuration when compared to the case 
without diquark condensation. For finite temperatures the masses are smaller 
than at T = 0. For asymptotically high temperatures and densities the EoS 
behaves like a relativistic ideal gas, where the relation pressure versus energy 
density is temperature independent. In contrast to the bag model where this 
behavior occurs immediately after the deconfinement transition, our model 
EoS has a temperature dependent P{e) relation also beyond this point. 

It has been shown for a hybrid star model which uses the quark matter EoS 
presented in this work that the possibility to obtain a stable star configuration 
with 2SC quark matter core depends on the formfactor of the quark interac- 
tion [42]. The Gaussian and Lorentzian formfactor models do allow a quark 
matter core, whereas the NJL formfactor model does not. Other NJL model 
calculations [26] confirm this conclusion. 

A mass defect of about 0.1 M & occurs in the cooling of a hot configuration 
to a cold one already without diquark condensation. It is mainly determined 
by the change of the configuration due to changes in the EoS. Cooling the 
star with diquark condensation results even in a lowering of the mass defect 
contrary to naive expectations. For a more realistic calculation of the cooling 
evolution one has to start from isentropic initial temperature profiles, which 
may change the size of the mass defect but not our main conclusions. 
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Table 1 

Parameters for the different formfactors models discussed in the text. Critical tem- 
peratures T c for the chiral phase transition at n q = and critical chemical potentials 
Hq(T = 0) for the onset of diquark condensation at T = 0. 
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40 


1.72 


9.67 
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1.553 
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8.44 


9.81 


1.893 


60 


1.61 


9.46 


7.95 


1.653 


1.45 


8.55 
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Table 2 

Characteristic quark star properties for Gaussian (left) and Lorentzian (right) form- 
factor models: maximal possible masses M max , radii R, central baryon densities 
ub{0) and central energy densities e(0) for different temperatures. 
*)Diquark gap is set to A = 0. 
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Fig. 1. Solutions for the chiral gap (p and the diquark gap A for different formfactors 
at T = in (3— equilibrium quark matter, G2 = G\. 
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Fig. 2. Solutions of the chiral gap <f> and the diquark gap A for the Gaussian model 
formfactor in the chiral limit (left panel) and for finite current quark mass tuq = 2.41 
MeV (right panel) at different temperatures T = 0, 40, 60 MeV for j3— equilibrium 
quark matter, G2 = G\. 
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Fig. 3. Phase diagrams for different formfactor modeis: Gaussian (solid lines), 
Lorentzian a = 2 (dashed lines) and NJL (dash-dotted) under conditions of 
/3-equilibrium and charge neutrality. The color-superconducting phase which ex- 
ists for G2 = G\ is absent for G2 = 0.75 G\. In the inset we show for the Gaus- 
sian model the comparison of the numerical result with the modified BCS formula 
T C A = 0.57 A(T = 0, fi q ) g(fiq) for the critical temperature of the superconducting 
phase transition. 
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Fig. 4. Dependence of the pressure on the baryochemical fis of quark matter for 
the Gaussian formfactor model for different temperatures T = 0, 40, 60 MeV. To 
demonstrate the effect of the stiffening of the EoS due to diquark condensation, we 
show the T = result for both A ^ (solid line, G2 = G\ ) and A = (dash-dotted 
line, G 2 = 0.75 G{). 
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Fig. 5. Matter distribution inside the quark star (lower panel) and accumulated mass 
(upper panel) for the quark stars with the same central baryon density n#(0) = 5 no 
and cut on ns{R) = 0.5 no, for different temperatures. The case A = is calculated 
for G 2 = 0.75 Gi, the others for G 2 = Gi. 
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Fig. 6. Stable configurations of quark stars for different temperatures T = 0, 40, 60 
MeV. Mass as a function of central baryon density (left panels) and of the radius 
(right panels) for the Gaussian model (upper graph) and for the Lorentzian model 
(lower graph). The case A = is calculated for G2 = 0.75 Gi, the others for 
G% = G\. 
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Fig. 7. Hot (T = 40 MeV) versus cold (T = 0) quark star configurations for the 
Gaussian model; case A: with diquark condensation (dashed versus full lines) and 
case B: without diquark condensation (dash-dash-dotted versus dash-dotted lines). 
When a quark star with initial mass Mj cools down from T = 40 MeV to T = at 
fixed baryon number Nb the mass defect AM occurs. The case A = is calculated 
for G 2 = 0.75 Gx, the others for G 2 = G X . 
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Fig. 8. Same as Fig. 7 for the Lorentzian model. 



24 



